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Abstract 



Polymerization dynamics of single actin filaments is investigated theoreti- 
cally using a stochastic model that takes into account the hydrolysis of ATP- 
actin subunits, the geometry of actin filament tips, the lateral interactions 
between the monomers as well as the processes at both ends of the polymer. 
Exact analytical expressions are obtained for a mean growth velocity and 
for dispersion in length fluctuations. It is found that the ATP hydrolysis 
has a strong effect on dynamic properties of single actin filaments. At high 
concentrations of free actin monomers the mean size of unhydrolyzed ATP- 
cap is very large, and the dynamics is governed by association/dissociation 
of ATP-actin subunits. However, at low concentrations the size of the cap 
becomes finite, and the dissociation of ADP-actin subunits makes a signif- 
icant contribution to overall dynamics. Actin filament length fluctuations 
reach the maximum at the boundary between two dynamic regimes, and 
this boundary is always larger than the critical concentration. Random 
and vectorial mechanisms of hydrolysis are compared, and it is found that 
they predict qualitatively similar dynamic properties. The possibility of at- 
tachment and detachment of oligomers is also discussed. Our theoretical 
approach is successfully applied to analyze the latest experiments on the 
growth and length fluctuations of individual actin filaments. 
Key words: ATP hydrolysis; actin filament growth 
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Introduction 



Actin filaments are major component of cytoskeleton in eukaryotic cells, and 
they play important roles in many biological processes, including the orga- 
nization of cell structures, transport of organelles and vesicles, cell motility, 
reproduction and endocytosis (jlj. Biological functions of actin fil- 

aments are mostly determined by the dynamic processes that take place 
during the growth or shrinking of these biopolymers. However, our under- 
standing of mechanisms of assembly and disassembly of these filaments is 
still very limited. 

In recent years the number of experimental investigations of the growth 



dynamics of rigid cytoskeleton filaments, such as actin filaments and mi- 
crotubules, at a single-molecule level have increased significantly (4, ^ ^, 



121 ). Dynamic behavior of individual microtubules have 
been characterized by a variety of experimental techniques such as video 
and electron microscopv, fluorescence spectroscopy, and optical trap spec- 



■on microscopv, nuc 

(IflflQJH), w 



trometry 
have just started (jlDl . 



11 



12 



lereas the studies of the single actin filaments 
. The assembly dynamics of individual actin 
filaments revealed a treadmilling phenomenon, i.e., the polymer molecule 
tends to grow at the barbed end and to depolymerize at the pointed end 



( 111). Similar picture has been observed earlier for microtubules 



<l], 



li). 



Although the conventional actin filaments do not exhibit the dynamic insta- 



bility as observed in microtubules 



it was shown recently that the 



DNA-segregating prokaryotic actin homolog ParM displays two phases of 



polymer elongation and shortening f|l^l. ATP hydrolysis is not required 
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for actin assembly (jTO), but it is known to play an important role in the 
actin polymerization dynamics. Experimental observations suggest that the 
nucleotide bound to the actin filament acts as a timer to control the filament 



turnover during the cell motility (|18l ) . Hydrolysis of ATP and release of in- 
organic phosphate are assumed to promote the dissociation of the filament 
branches and the disassembly of ADP-actin filaments ( 2) • 



1^ 



Recent experimental studies of the single actin filament growth 
have revealed unexpected properties of actin polymerization dynamics. A 
large discrepancy in the kinetic rate constants for actin assembly estimated 
by average length change in the initial polymerization phase and determined 
from the analysis of length fluctuations in the steady-state phase (by a factor 



of 40) has been observed. Several 



servation has been proposed 



m)ssib 



12 



3, 



e explanations of this intriguing ob- 



211). First, the actin polymerization 



dynamics might involve the assembly and disassembly of large oligomeric 
actin subunits. However, this point of view contradicts the widely accepted 



2lh . In addition, as we 



picture of single-monomer polymerization kinetics (|l 
argue below, it would require the association/dissociation of actin oligomers 
with 30-40 monomers, but the annealing of such large segments has not 



been observed in experiments or has been excluded from analysis (jll 



1^. 



Second, stochastic pauses due to filament-surface attachments could increase 
the apparent dispersion in length of single actin-filaments, although it seems 



that the effect is not significant 



Third, the errors in the experimental 



measurements could contribute into the observation of large apparent diffu- 
sion constants (|l^ . Another possible reason for the discrepancy is the use 
of oversimplified theoretical model in the analysis that neglects the polymer 
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structure and the lateral interaction in the actin filament. However, the 
detailed theoretical investigation of the growth of single actin filaments 
indicates that large length fluctuations still cannot be explained by correctly 
describing the structure of the filament's tip and the lateral interactions be- 
tween the monomers. 

The fact that hydrolysis of ATP bound to the actin monomer is im- 
portant stimulated a different model to describe the actin polymerization 
dynamics (21). According to this approach, the ATP-actin monomer in 
the filament can be irreversibly hydrolyzed and transformed into the ADP- 
actin subunit. The polymer growth is a process of adding single ATP-actin 
monomers and deleting hydrolyzed or unhydrolyzed subunits, and the actin 
filament consists of two parts - a hydrolyzed core in the middle and unhy- 
drolyzed caps at the polymer ends. Large length fluctuations are predicted 
near the critical concentration for the barbed end. Although this model 
provides a reasonable description of the filament growth rates for different 
ATP-actin concentrations, the position of the peak in dispersion is below the 



critical concentration for the barbed end, while in the experiments , 



1^ 



large dispersion is observed at, or slightly above, the critical concentration. 
In addition, the proposed analytical method (without approximations) (2lll 
cannot calculate analytically cap sizes and dispersion. 

The goal of this work is to develop a theoretical model of polymeriza- 
tion of single actin filaments that incorporates the ATP hydrolysis in the 
polymer, the structure of the filament tips, lateral interactions between the 
monomers and dynamics at both ends. Our theoretical method is based on 
the stochastic models developed for describing the growth dynamics of rigid 
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multifilament biopolymers 



221 ). and it allows to calculate explicitly all 
dynamic parameters of the actin filament's growth. Different mechanisms of 
ATP hydrolysis in actin filaments are compared. The possibility of adding or 

yz e the 



1^. 



deleting oligomeric actin subunits is also discussed. Finally, we anal' 
latest experiments on the growth dynamics of single actin filaments (jlll . 

This article is organized as follows. The model of the actin polymer- 
ization is presented in the Section II. The effect of different mechanisms of 
hydrolysis on the actin filament growth is discussed in the Section III. The 
possibility of assembly and disassembly of oligomers is critically evaluated 
in the Section IV. The summary and conclusions are given in the Section V, 
while the technical details of calculations can be found in the Appendix. 

Model of actin filament assembly 

Let us consider an actin filament as a two-stranded polymer, as shown in 
Fig. 1. It consists of two linear protofilaments. The size of the monomer 
subunit in this polymer is equal to d = 5.4 nm, and two protofilaments are 
shifted with respect to each other by a distance a = d/2 = 2.7 nm (1, 
Each monomer in the actin filament lattice carries a nucleotide molecule: 
it can be either ATP, or ADP (see Fig. 1). Shortly after actin monomers 
assemble into filaments, the ATP is hydrolyzed to ADP. For simplicity, we 
neglect the intermediate states of hydrolysis for actin monomers when both 
products of hydrolysis, ADP and Pj (inorganic phosphate), are bound to the 
actin monomer. Thus, we only consider two states of actin monomers in the 
filament - hydrolyzed and unhydrolyzed. It is argued that only the dynam- 
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ics of capped (unhydrolyzed) and uncapped (hydrolyzed) states effect the 

n 

length fluctuations in the actin filaments (21). The hydrolyzed nucleotide 
remains bound to the polymer and at the physiological conditions it is not 
exchangeable with free ATP molecules from the solution. The dynamical 
and biochemical properties of ATP-bound (T-state) and ADP-bound (D- 
state) monomers are known to be different ([2^. The dissociation rate of 
ADP-actin subunits from the actin filaments is estimated to be 2-5 times 
larger than the rate for the ATP-actin subunits, whereas the association rate 
is considerably slower (by a factor of 10) than that for the non-hydrolyzed 



analog (jl2 



23). 



ATP hydrolysis plays an important role for the overall actin filament 
assembly dynamics, however the details of this process are not clear. Sev- 
eral mechanisms of ATP-actin hydrolysis in the filament have been pro- 

nriQ 

posed. In a random mechanism ()24l . 1251 . \2ai ) any ATP-actin subunit can 
hydrolyze in a stochastic manner independently of the states of the neigh- 
boring monomers. The rate of hydrolysis in this case is proportional to the 
amount of non-transformed nucleotide in the polymer. A different approach 
is a sequential, or vectorial, mechanism ( 27 . 12 J [2^, that assumes a high 
degree of cooperation during the hydrolysis. According to this mechanism, 
the recently assembled actin monomer hydrolyzes its ATP only if it touches 
the more interior already hydrolyzed subunits. In this mechanism, there is 
a sharp boundary between the unhydrolyzed cap and the hydrolyzed core 
of the filament, while in the random mechanisms there are many interfaces 
between ATP-actin and ADP-actin subunits. Finally, it is also possible that 
a "mixed" mechanism, that combines the properties of random and vecto- 
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rial approaches, describes the hydrolysis in actin filaments. The available 
experimental data cannot clearly distinguish between these mechanisms. In 
our model we assume a vectorial mechanism, although, as we show below, 
the exact details of hydrolysis do not influence much the dynamic properties 
of the actin filament's growth. 

There are infinite number of possible polymer configurations depending 
on the nucleotide state of each monomer and the geometry of polymer ends 
((2^: see Fig. 1. However, we assume that only the so-called "one-layer" 
configurations, where the distance between two edge monomers at parallel 

ymerization dynam- 



22I ). that showed 



protofilaments is less than d, are relevant for actin po 
ics. This is based on the previous theoretical studies (|2nl . 
that the one-layer approach is an excellent approximation to a full dynamic 
description of growth of two-stranded polymers with large lateral interac- 
tions between the subunits. It is also known that for actin filaments the 
lateral interaction energy is larger than 5 k^T per monomer (2^, ^), and 
it strongly supports the "one-layer" approximation. 

Each configuration we label with two pairs of integers, (/i, fci; ^2); 
where li is the total number of monomers (hydrolyzed and unhydrolyzed) in 
the i-th protofilament, while ki specifies the number of ATP-actin subunits 
in the same protofilament. For example, the configuration A from Fig.l is 
labeled as (2,1;3,2), while the configuration B is described as (2,1;3,1). The 
polymerization dynamics at both ends of the actin filament is considered 
independently from each other. 

As shown in Fig. 1, at each end free ATP-actin molecules from the 
solution can attach to the actin filament with the rate u = k^c, where kx 
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is the ATP-actin polymerization rate constant and c is the concentration 
of free ATP-actin species in the solution. Because of the excess of free 
ATP molecules in the solution only ATP-actin monomers are added to the 



filament ((27 



28, 3). 



We also assume that the dissociation rates of actin 
monomers depend on their nucleotide state, and only the leading subunits 
dissociate from the filament. Specifically, ATP-actin monomer may detach 
with the rate wt, while the hydrolyzed subunit dissociates with the rate wd- 
see Fig. 1. In addition, the sequential vectorial mechanism of hydrolysis is 
assumed, i.e., ATP-actin monomer can transform into ADP-actin state with 
the rate if it touches two already hydrolyzed subunits. 

The growth dynamics of single actin filaments can be determined by 
solving a set of Master equations for all possible polymer configurations. 
The mathematical derivations and all details of calculations are given in 
Appendix. Here we only present the explicit expressions for the dynamics 
properties of actin filament growth at stationary state. Specifically, the 
mean growth velocity is equal to 



V = ^{u - WTQ - wnil - q)}, (1) 



and dispersion is given by 



J 2{wd-wt){u + wdc 
D = — <u + wxq + wd{1 - qj -\ 



WT + rh 
for < g < 1, where 
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The parameter q plays a critical role for understanding mechanisms of 
actin growth dynamics. It has a meaning of probability that the system is 
in a "capped" state with Neap > 1 ATP-actin monomers , i.e., it is a fraction 
of time that the actin filament can be found in any configuration with at 
least one unhydrolyzed subunit. For example, in Fig. 1 the configurations 
A and B are capped (with Neap = 3 and 2, correspondingly), while the 
configuration C is uncapped with Neap = 0. The parameter q increases 
linearly with the concentration of free ATP-actin monomers because of the 
relation u = Utc. However, it cannot be larger than 1, and this observation 
leads to existence of a special transition point with the concentration 



Above the transition point we have q{c > c') = 1, and the probability to 
have a polymer configuration with Neap = is zero and the unhydrolyzed 
ATP cap grows steadily with time. At large times, for c > c' the average 
length of ATP cap is essentially infinite, while below the transition point 
(c < c') this length is always finite. 

At the critical concentration for each end of the filament, by definition, 
the mean growth velocity for this end vanishes. Using Eqs. Q and El it can 
be shown that 



An important observation is the fact the critical concentration is always 



(4) 



wp {WT + rh) 
kr {wD + rh) ' 



(5) 
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below the transition point, 

CcTit = \ c . (6) 

wd + rh 

Because at concentrations larger than the transition point the dissocia- 
tion events of hydrolyzed actin monomers are absent, the explicit expressions 
for the mean growth velocity and dispersion in this case are given by 



The calculated mean growth velocity for the barbed end of the actin filament 
is shown in Fig. 2 for parameters specified in Table 1. It can be seen that the 
velocity depends linearly on the concentration of free ATP-actin particles 
in the solution, although the slope changes at the transition point. This 
is in agreement with experimental observations on actin filament's growth 
However, the behavior of dispersion is very different — see Fig. 3. 
It also grows linearly with concentration in both regimes, but there is a 
discontinuity in dispersion at the transition point. From Eqs. and (O 
we obtain that the size of the jump is equal to 



The origin of this phenomenon is the fact that ADP-actin subunits disso- 
ciations contribute to overall growth dynamics only below the transition 
point c'. Note, as shown in Fig. 3, this contribution can increase the length 
fluctuations when wd > wt (the barbed end of the filament), or dispersion 
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TABLE 1 Summary of rate constants 

Rate constant Reaction Barbed end Pointed end Source 

kr, IJ-M-^ s"^ ATP-actin association IL6 L3 Refs. (18, 23) 

WT,s-^ ATP-actin dissociation 1.4 0.8 Refs. (18, 23) 

WD,s-^ ADP-actin dissociation 7.2 0.27 Refs. (18, 23) 

rh, s'l ATP hydrolysis 0.3 0.3 Ref. (26) 



can be reduced for wd < wt (the pointed end of the filament). The jump 
disappears when wd = wt- For the barbed end of the actin filament we cal- 
culate, using the parameters from the Table 1, that D{c')/ Dq{c') ~ 20.6 and 
it approaches to 26.5 when rh — s- 0. This result agrees quite well with the 
experimentally observed "apparent" difference in the kinetic rate constants 



(35-40 times) (ju, 



1^. 



In order to compare our theoretical predictions with experimental obser- 
vations the dynamics at both ends should be accounted for. However, as we 
showed earlier the total velocity of growth and the overall dispersion 
are the sums of the corresponding contributions for each end of the filament. 
The parameters we use in the calculations are shown in Table 1. Experi- 
mental measurements of actin filament's growth suggest that the hydrolysis 
rate is rather small, and we took = 0.3 s"^ as given in the most recent 
investigation Note, however, that this value is determined assuming 

the random mechanism of hydrolysis. But the specific value of r^ does not 
strongly influence our calculations. More controversial is the value of ADP- 
actin dissociation rate constant wd- Most experiments indicate that wd is 
relatively large, ranging from 4.3 s~^ to 11.5 s~ ^ however the lat- 
est measurements performed using FRET method (l^) estimated that the 
dissociation rate is lower, wd = 1-3 s^^. We choose for wd the value of 7.2 
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as better describing the majority of experimental work. 

For the actin filament system with the parameters given in Table 1 we 
can calculate from Eq. that the critical concentration for the barbed end 
is Ccrit — 0.141 /^M, while for the pointed end it is equal to Ccrit 

~ 0.401 

/iM. However, the contribution of the pointed end processes to the overall 
growth dynamics is very small. As a result, the treadmilling concentration, 
when the overall growth rate vanishes, is estimated as ctm — 0.144 /xM, 
and it is only slightly above the critical concentration for the barbed end 
(see Fig. 3). The treadmilling concentration also almost coincides with the 
transition point for the barbed end, as can be calculated from Eq. ^ c' ~ 
0.147 /xM. According to Eq. 121 the dispersion at treadmilling concentration 
at stationary-state conditions is equal to D{ctm) — 31.6 sub'^-s^^, with the 
contribution from the pointed end equal 0.5%. From experiments, the values 
29 sn62.s-i (ll and 31 sub'^-s'^ (d) are reported for the filaments grown 
from Mg-ATP-actin monomers, and 25 sub'^-s^^ (11) is the dispersion for 
Ca-ATP-actin filaments. The agreement between theoretical predictions 
and experimental values is very good. It is also important to note that, 
in contrast to the previous theoretical description (21), our model predicts 
large length fluctuations slightly above the Ccrit for the barbed end of the 



1^ 



filaments, exactly as was observed in the experiments (jlll . 

The presented theoretical model allows to calculate explicitly not only 
the dynamic properties of actin growth but also the nucleotide composition 
of the filaments. As shown in Appendix, the mean size of the cap of ATP- 
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actin monomers is given by 

< N,ap >= = -J^. (9) 

I — q c— c 

Then at the critical concentration for the barbed end, Ccrit 

~ 0.141 fiM, 

the cap size at the barbed end is N^ap — 24, while the cap size at the 
pointed end at this concentration (with the transition point c' ~ 0.846 ^M) 
is less than 1 monomer. These results agree with the Monte Carlo computer 
simulations of actin polymerization dynamics (21). Large ATP cap appears 
at the barbed end of actin filament and smaller cap is found at the pointed 
end This is quite reasonable since the transition point for the barbed 

end is much smaller than the corresponding one for the pointed end. The 
overall dependence of Neap on the concentration of free actin monomers is 
shown in Fig. 4. 



Random vs. vectorial ATP hydrolysis in actin fila- 
ments. 

An important issue for understanding the actin polymerization dynamics 
is the nature of ATP hydrolysis mechanism. In our theoretical model the 
vectorial mechanism is utilized. In order to understand what features of 
actin dynamics are independent of the details of hydrolysis, it is necessary 
to compare the random and the vectorial mechanisms for this process. In our 
model of the polymer's growth with vectorial mechanism, the actin filament 
consists of two parallel linear chains shifted by the distance a = d/2 from 
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each other. According to our dynamic rules, the addition (removal) of one 
actin subunit increases (decreases) the overall length of the filament by the 
distance a. Then the growth dynamics of two-stranded polymers can be 
effectively mapped into the polymerization of single-stranded chains with 
an effective monomer's size dgff = d/2. 

A single-stranded model of the actin filament's growth that assumes 
association of ATP-actin monomers and dissociation of ATP-actin and ADP- 
actin subunits along with the random hydrolysis has been developed earlier 



meaning of the probability to find the leading subunit of the polymer in 
the unhydrolyzed state. However, the parameter q in the random hydrolysis 
model has a more complicated dependence on the concentration than in the 
vectorial model. It can be found as a root of the cubic equation, 

[u- {u + WT + rh)q] {u - wtq) + [wtQ + wd(1 - q)] [u - {wt + rh)q] g' = 0, 



with the obvious restriction that < (7 < I. 

Using the parameters given in the Table 1, the fraction of the capped 
configurations for two different mechanisms of hydrolysis is shown in Fig. 
5. The predictions for both mechanisms are close at very low and very 
high concentrations, but deviate near the critical concentration. It can be 
understood if we analyze the special case wt = wd, although all arguments 
are still valid in the general case of wd 7^ wt- From Eq. ^Jwe obtain 




In this model the parameter q is also introduced, and it has a 



(10) 




(11) 
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In the limit of very low concentrations the fraction q approaches 

u 

WT + rh 

while for c ^ 1 it can be described as 



1 ^ . (13 

u 



Generally, in the limit of very low hydrolysis rates the random and the 
vectorial mechanisms should predict the same dynamics, as expected. It 
can be seen by taking the limit of r/j — > in Eq. Ilfl( which gives q = u/wt 
for u < Wt and q = 1 for u > wt- These results are illustrated in Fig. 5. 

In order to calculate the size of the ATP-cap in the actin filament for 
the random mechanism we introduce a function P„ defined as a probability 
to find in the ATP-state the monomer positioned n subunits away from the 
leading one. Then it can be shown that this probability is exponentially 



decreasing function of n 



%i = l-^:^, P,=q. (14) 
Fn u- wrq 



The size of the unhydrolyzed cap in the polymer is associated with the total 



number of ATP-actin monomers (|21l). 

oo 

The results of the different mechanisms for N^ap are plotted in Fig. 4. The 
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random and vectorial mechanism agree at low concentrations, but the pre- 
dictions differ for large concentrations. Only in the limit of small hydrolysis 
rates the predictions from the two hydrolysis mechanisms start to converge. 

The mean growth velocity in the model with the random mechanism is 
given by 

V = deff[u - WTQ - wd{1 - q)] , (16) 

which is exactly the same as in the vectorial mechanism (see Eq. ^ , although 
the fraction of ATP-cap configurations q have a different behavior in two 
models. Mean growth velocities for different mechanisms are compared in 
Fig. 2. Again, the predictions for different mechanisms of hydrolysis con- 
verge for very low and very high concentrations of free actin, but differences 
arise near the critical concentration. 

In the model with the random mechanism of hydrolysis (0, Q the 



analytical expressions for dispersion have not been found. However, from 
the comparison of the fraction of ATP-capped configurations q, the size of 
the ATP-cap Neap and the mean growth velocity V it can be concluded that 
both mechanisms predict qualitatively and quantitatively similar picture for 
the dynamic behavior of the single actin filaments. Thus, it can be expected 
that, similarly to the vectorial mechanism, there is a peak in dispersion 
near the critical concentration in the model with the random hydroly sis, in 
agreement with the latest Monte Carlo computer simulations results (|2l|). 
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Assembly /disassembly of oligomers in actin filament 
dynamics 

The association and dissociation of large oligomers of actin monomers has 
been suggested as a possible reason for large fluctuations during the elonga- 
tion of single actin filaments ((ill . Il2t ) . Let us consider this possibility more 
carefully. Suppose that the oligomeric particles that contain n ATP-actin 
monomers can attach to or detach from the filament. Then the mean growth 
velocity can be written as 



Ti n 

V{n) = —[u{n)-WT{n)l (17) 



where u{n) and WT{n) are the assembly and disassembly rates of oligomeric 
subunits. Similarly, the expression for dispersion is given by 

D{n) = ^—^ [u{n) + WT{n)] . (18) 



At the same time, in the analysis of the experimental data (|lll . Il2t ) the 
addition or removal of single subunits has been assumed. It means that the 
rates has been measured using the following expression, 

V = ^[u'ff -w'J^]. (19) 

Comparing this equation with Eq. I17[ it yields the relation between the 
effective rates u*^-^-^ and Wjf-^ per monomer and the actual rates n(n) and 
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wxiji) per oligomer, 

yfiff — nu{n), = nwT{n). (20) 

The substitution of these effective rates into the expression for dispersion 
(|18j) with n = 1 produces 

D{n = l) = — u'^f^ + Wjf^ = [u{n)+WT{n)] = ^^. (21) 
8 L J 8 n 

It means that dispersion calculated assuming the association/dissociation of 
monomers underestimates the "real" dispersion in n times, but not in 



suggest that only the addition or dis- 



The experimental results 
sociation of oligomers with n = 35 — 40 can explain the large length fluctua- 
tion in the single actin filaments if one accepts the association/dissociation 
of oligomers. These particles are quite large by size (~100 nm), and, if 
present in the system, they would be easily observed in the experiments. 
However, no detectable amounts of large oligomers have been found in stud- 
ies of kinetics of the actin polymerization. It has been reported (js^) that 
only small oligomers (up to n = 4 — 8) may coexist with the monomers and 
the polymerized actin under the special (not physiological) solution condi- 
tions. Therefore, it is very unlikely, that the presence of very small (if any) 
amounts of such oligomers might influence the dynamics and large length 
fluctuations in the actin growth. 
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Summary and conclusions 

The growth dynamics of single actin filaments is investigated theoretically 
using the stochastic model that takes into account the dynamics at both ends 
of filament, the structure of the polymer's tip, lateral interactions between 
the protofilaments, the hydrolysis of ATP bounded to the actin subunit, and 
assembly and disassembly of hydrolyzed and unhydrolyzed actin monomers. 
It is assumed that sequential (vectorial) mechanism of hydrolysis controls 
the transformation of ATP-actin subunits. Using the analytical approach, 
exact expressions for the mean growth velocity, the length dispersion, and 
the mean size of fluctuating ATP-cap are obtained in terms of the kinetic 
rate constants that describe the assembly and disassembly events, and the 
hydrolysis of nucleotides. It is shown that there are two regimes of single 
actin filament's growth. At high concentrations the size of the ATP-cap is 
very large and the the fully hydrolyzed core is never exposed at filament's 
tip. As a result, the disassembly of ADP-actin subunits does not contribute 
to the overall dynamics. The situation is different at low concentrations, 
where the size of ATP-cap is always finite. Here the dissociation of both 
hydrolyzed and unhydrolyzed actin monomers is critical for the growth dy- 
namics of filaments. The boundary between two regimes is defined by the 
transition point, that depends on the association/dissociation rate constants 
for the ATP-actin monomers and on the hydrolysis rate. For any non-zero 
rate of hydrolysis the transition point is always above the critical concen- 
tration where the mean growth velocity becomes equal to zero. The most 
remarkable result of our theoretical analysis is the non-monotonous behavior 
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of dispersion as the function of concentration and large length fluctuations 
near the critical concentration. These large fluctuations are explained by 
alternation between the relatively slow assembly /disassembly of ATP-actin 
subunits and the rapid dissociation of ADP-actin monomers. 

For the experimentally determined kinetic rates our theoretical analysis 
suggests that the contribution of the dynamics at the pointed end of the actin 
filament is very small. The treadmilling concentration for the system, as well 
as the transition point for the barbed end, are only slightly above the critical 
concentration for the barbed end. At this transition point the dispersion of 
length reaches a maximal value, and it is smaller for larger concentrations. 
Our theoretical predictions are in excellent agreement with all available ex- 
perimental observations on the dynamics of single actin filaments. However, 
more measurements of dispersion and other dynamic properties at different 
monomeric concentrations are needed in order to check the validity of the 
presented theoretical method. 

Since the exact mechanism of ATP hydrolysis in actin is not known, 
we discussed and compared the random and the vectorial mechanisms for 
the simplified effective single-stranded model of actin growth. It was shown 
that the mean growth velocity, the fraction of the configurations with unhy- 
drolyzed cap and the size of ATP-cap are qualitatively and quantitatively 
similar at low and at high concentrations, although there are deviations near 
the transition point. It was suggested then that the length fluctuations in 
the random mechanism, like in the vectorial mechanism, might also exhibit 
a peak near the critical concentration. Future experimental measurements 
of dynamic properties at different concentrations might help to distinguish 
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between the hydrolysis mechanisms in the actin growth. 

The possibihty of attachment and detachment of large oligomers in the 
single actin filaments has been also discussed. It was argued that the ex- 
perimental observations of large length fluctuations can only be explained 
by addition of oligomers consisting of 35-40 monomers. However, these 
oligomers have a large size and such events have not been observed or have 
been excluded from the analysis of the experiments. Thus the effect of 
the oligomer assembly /disassembly on the single actin filament's growth is 
probably negligible. 

Although the effect of ATP hydrolysis on polymerization dynamics of 
actin filaments has been studied studied before to best of our 

knowledge, the present work is the first that provides rigorous calculations 
of the mean growth velocity, dispersion, the size of ATP-cap and the fraction 
of capped configurations simultaneously. It is reasonable to suggest that this 
method might be used to investigate the dynamic instability in microtubules 
because the polymer can be viewed as growing in two dynamic phases. In one 
phase the ATP-cap is always present at the end of the filament, while in the 
second phase it is absent. Similar approach to investigate the dynamic phase 
changes has been proposed earlier |3g). The model can also be improved by 
considering the intermediate states of hydrolysis and the release of inorganic 
phosphate (3, 21), and the possible exchange of nucleotide at the terminal 
subunit of the barbed end of the actin filaments 
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Appendix: One-layer polymerization model with 
sequential ATP hydrolysis for two-stranded poly- 
mers. 

Let us define a function P{li,ki;l2,k2;t) as the probability of finding the 
two-stranded polymer in the configuration (/i, ki;l2,k2). Here li, ki = 0,1, ... 
(ki = ! or 2) are two independent parameters that count the total 

number of subunits (/j) and the number of unhydrolyzed subunits (ki) in 
the i-th protofilament. We assume that the polymerization and hydrolysis 



22). It 



in the actin filament can be described by "one-layer" approach (|2C . 
means that I2 = h or I2 = h + 1, and k2 = ki or k2 = ki ± 1 (see Fig. 1). 
Then the probabilities can be described by a set of master equations. For 
configurations with li = I2 = I and 1 < A; < / we have 

dP{l,k^l,k;t) ^ ^^^^ _ 1^ ^ _ 1. ^. + ^^p^i^ k;l + l,k + l;t) 

+rhP{l, k + l;l, k; t)-{u + WT + rh)P{l, k; I, k; t), (Al) 
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and 

dP{l,k;l,k - l;t) 



= uP{l -l,k-l;l,k-l-t) + wtP{1, k;l + l, k; t) 



dt 

+rhP{l, k; I, k; t)-{u + WT + rh)P{l, k;l,k- 1; t). (A2) 

Similarly for the configurations with li = I2 — 1 = I and 1 < A; < / + 1 the 
master equations are 

dPi.hk-l^l + lMt) ^ ^p^^^ j,_i.i^k-l-t) + wtP{1 + l,k-l + l, k; t) 

+rhP{l, k-l + 1, k; t)-{u + WT + rh)P{l, k-l;l + l,k; t)iA3) 

and 

^^^^^^^i^^^^ = uP(Z, fc; Z, A; - 1; + w;tP(Z + 1, )fc + 1; / + 1, it; t) 

+rhP{l, k;l + l,k + l;t) - {u + WT + rh)P{l, k;l + l, k; t). (A4) 

Then the polymer configurations without ATP-actin monomers {k = 0) can 
be described by 

dP(l,0;l,0;t) ^ , . , s , ^ x 

— ^ ' ' ' ^ = wrPil, 0; Z + 1, 1; t) + wdP{1, 0; / + 1, 0; t) 

-{u + WD)P{l,0;l,0;t), (A5) 

and 

"^^^^'"'ir^'"'^^ " ""^^^^ + 1, 1; / + 1, 0; t) + wdP(.1 + 1, 0; Z + 1, 0; t) 

-{u + WD)P{l,0;l + l,0;t). (A6) 
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Finally, for the configurations consisting of only unhydrolyzed subunits we 
have 



dP{l,l;l,l;t) 
di 



uP{l I; t) + wtP{1, l;l + 1,1 + l;t) 

-{u + WT + rh)P{l,l;l,l;t), (A7) 



and 



dP{l,l;l + 1,1 + l;t) 
di 



uP{l, I; I, I; t) + wtP{1 + 1,1 + l;l + 1,1 + l;t) 

-{u + WT + rh)P{l, + + t), (A8) 



The conservation of probability leads to 



E 



+ 00 



^ P{1, k; I, k;t)+^ P{1, k + l;l, k; t) 



1=0 



1=0 



+ J2 ^(^ A:; / + 1, A; + 1; t) + ^ P{1, k;l + l, k; t) 



1=0 



1, 



(A9) 



at all times. 



0), 



Following the method of Derrida (!40'), we define two sets of auxiliary 
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functions {k = 0, 1, ...), 

+ 00 

l^k 

+ 00 

1=0 

Bl^+^{t) = Y,P{W + l,k + l-t), 

l^k 

1=0 

+00 

SO(t) = ^P(Z,/;Z,/;i), 

1=0 

^^(t) = ^P(Z,/;Z + l,Z + l;i); (AlO) 

1=0 

and 



+ 00 ^ 

Cl,{t) = Y,il + -)P{l,k;l,k;t), 

+00 1 

1=0 

+00 

Cfe,ik+i(i) = E + k;l + l,k + l;t), 

+ 00 

ci(t) = E + ^; ^ + 1> ^ + 1; (All) 
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Note that the conservation of probability gives us 
+00 

Kkit) + Bl+^.k{t) + + Blkit) + B\t) + B\t) = 1. (A12) 

k=0 

Then from the master equations ()A1|) . ()A2|) . ()A,3|) and ()A4|) we derive for 
k > 1 

= nBU,kii) + ^TBlkA^) + rnBl^, ,{t) -{U + WT + rH)Bl,{t), 

^^^^^ = ^B^k-i,k-iit) + WTBl^it) + rnBl^it) - {u + wt + rh)Bl_, ,{t), 
dBl.,{t) 



dt 



^Bl_,.,{t) + wrBUt) + rhBlJt) - (u + wt + rh)Bl,_,{t), 



dB^ (t) 

= ^Bl,_,{t) + WTBl^,^,{t) + rj,Bl,^,it) -iu + WT + r,)i?i,,((AX3) 



while the master equations (|A5|1 . (|A6|) for k = yield 
dBUt) 



dt 

dBUt) 



wrBUt) + woBUt) + rhB%{t) -{u + WD)BUt), 



TB%it) + WDB°^o{t) + ThBl^{t) -{u + WD)Bloit)i^'^^) 



dt 

Finally, equations HA7|) and HA8|) lead to 
dB^{t) _ 

{u + WT)B^{t) -{U + WT + rh)Bi{t), (A15) 



- iu + WT)B^t) - {u + WT + rh)B^{t), 
dt 

dB^{t) 



dt 

Similar arguments can be used to describe the functions C^^, C^j^^^., 
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^i-fe+i' ^kk- Specifically, for > 1 we obtain 

= ^C'Ukit) + ^TClk+iit) + rkC'k+i,kit) -{n + WT + r,)Cl,{t) 
m^Bl_,^kit)-^TBlk+iit)h (A16) 

dC^ (t) 

= ^C'k-i,k-i{t) + ^TCl,{t) + rf,Cl,{t) -{U + WT + rh)Cl_,^,{t) 
H[uBl_,,_,{t)-WTBl,{t)], (A17) 



dC^ (t) 

= ^Cl_,^k-iit) + ^TClkit) + rhCl,{t) -{U + WT + ru)Cl,_,{t) 

muBl_,,_,it)-WTBlJt)], (A18) 



= ^Cl,_,{t) + WTC^,^,^,{t) + rnCl,^,{t) -{u + wt + r,)Cl,{t) 
m^Bl,_,{t)-WTB^,^,^,{t)]. (A19) 

For k = the expressions are the following, 
dC^ (t) 

— ^ = wrCl.it) + woCl^it) + mCUt) -iu + wd)C^oA^) 

-^[wTBl,it)+WDBl,it)], (A20) 



dC^ (t) 

— ^ = wrCloit) + WDC%{t) + rnCl,{t) - (n + WD)Cloit) 

_ 1 U,,„ rO (+\ I rO 



^[wTB%{t) + WDBlo{t)]. (A21) 
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Again following the Derrida's approach we introduce an ansatz that 
should be valid at large times t, namely, 

BlrrXt) - blm, ^ + U^m = 0, 1; \k - m\ < 1). (A22) 

At steady state dBi^^{t)/dt = 0, and Eas. 1X131 and IXTl yield for A: > 1 

= + ^Tblk+i + ^hbl+^^k -{U + WT + rh)bl i^, 

= + WTb^k + ^hblk -{n + WT + rh)bl_^^k^ 

= + '^Tb\^k + ^hblk -{U + WT + rh)blk-v 

= ^bl,_, + WTbl^,^„ + r,,bl,^, -iu + WT + r06^,fc, (A23) 

while for = we obtain 

= wrbl i + wnbl o + rhb^ Q - {u + WD)bQ Q, 

= WTb% + WDb^o + rhbli - {u + wdH^q. (A24) 

Finally, from Eq. IA15I we have 

= {u + WT)b^ -{u + WT + rh)b°, 

= (^u + WT)b° - {u + WT + rh)bi. (A25) 

Due to the symmetry of the system we can conclude that the probabilities 
bkk~bkk ^^'^ ^fe+i k ~ bk k+i- Then the solutions of Eas. IA231 and IA2I1 can 
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be written in the following form, 

bl+,,u = biu+i = W-^)<f'^\ (A26) 

where A: = 0, 1, .., and 

q = < 1. (A27) 

WT + Th 

In addition, Eqs. lA25l have only a trivial solution = = 0. Recall that b^ 
and b^ give the stationary-state probabilities of the polymer configurations 
with all subunits in ATP state, i.e., it corresponds to the case of very large 
k. The solution agrees with the results for b^ ^. and 6^ at A; — > oo (see 

For q > 1 the systems of equations (|A23|) and HA24|) have the trivial 
solutions, with all b\^ = Q for finite k and m. It means that at the stationary 
conditions the polymer can only exist in the configurations with very large 
number of unhydrolyzed subunits and the size of ATP-cap is infinite, while 
for g < 1 the size of ATP-cap is always finite. The case g = 1 is a boundary 
between two regimes. At this condition there is a qualitative change in the 
dynamic properties of the system. 

To determine the coefficients a\ ^ and ^ from Eq. IA221 the ansatz 
for the functions ^ is substituted into the asymptotic expressions 
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IA2ip . yielding for /c > 1, 



= u4-i,k + ^Talk+i + rhal+^,k -{U + WT + rh)al k, 
= ^a°-i,fc-i + ^Tal^k + ^hal j, -{u + wt + rh)al_^^j,, 
= W-i,fe-i + wral^k + ^hO^Kk -{u + wt + rh)al^k-i^ 
= + "^Tffl^+i.fc + r^a^^fc+i - (n + 'u;t + r,,)ai , (A28) 

At the same time, for A; = we obtain 

= WTO? + ^t'Dao,o + ^/iOq,! - + WD)alfi. (A29) 
The coefficients satisfy the following equations (for /c > 1), 



+h\<-i,k-WTb\k+il (A30) 

= uT^-i,k-i + ujtTI^^ + ThTlt, -{U + WT + rh)Tl_-^ ^^ 

+\[<-i,k-i-^Tblki (A31) 

mubl_,,_,-WTblA (A32) 
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«U = uTl^k_^ + wtT^_^^^^ + rhTl j^_^-^ -{U + WT + rh)Tlk 

+l[<,k-i-^Thl^i,kl (A33) 

For A: = the expressions are given by 

ao,o = wtTI^ + wdTI^ + ThTl^ - (n + wd)TIq 

-\[wThl^ + WDblol (A34) 



00,0 = "f^TTi^o + woT^fi + rhTl^ - (n + wd)?}} o 
-\[wTh% + WDhl^Q]. 



(A35) 



Comparing Eas. IA231 and 1^241 with expressions IA28I and IA291 we con- 
clude that 



^\,m — ^^fc,m' — 0, 1), 



(A36) 



with the constant A. This constant can be calculated by summing over the 
left and right sides in Eq. IA36I and recalling the normalization condition 
(|A12|) . The summation over all at^i in Eqs. IA281 and 1X291 produces 



fc=0 



7° 



\ [{u - wt) - (wD - 'wt)(6o,o + ^o,o)] 



(A37) 



To determine the coefficients TJ^^ , we need to solve Eqs. rO?M35l 
Again, due to the symmetry, we have T^^ = T^^^ = T2k, and T^_^_l f. = 
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= T2k+i for all k. The solutions for these equations are given by 

n = + (A38) 

4 WT + Th 



where A; = 0, 1, ... and Tq is an arbitrary constant. 

It is now possible to calculate explicitly the mean growth velocity, V, 
and dispersion, D, at steady-state conditions. The average length of the 
polymer is given by 

( +00 +00 +00 +00 

k=0 1=0 k=0 1=0 

+00 +00 +00 +00 

+ E Eil + '^)Pil,k;l,k + l;t)+Z Y.{l + ^)P{W + lMt] 

k=0 1=0 k=0 1=0 

= dE Cl,{t) + C',^,^,{t) + Cl,_^,{t) + Cl,{t) . (A39) 

k=0 '- 

Then, using Eq. IA361 we obtain for the velocity 



hoo 



V = /nn - < l{t) >= + bk+i,k + Yl ^U+i + E ^ik ) = dA. 

\fe=0 k=0 k=0 k=0 / 

(A40) 



A similar approach can be used to derive the expression for dispersion. 
We start from 



( +00 +00 +00 +00 

E E a + h)^Pil, k; /, fc; t) + E E a + l?P{l, k + l;l, k; t) 
k=0 1=0 k=0 1=0 

+00 +00 +00 +00 \ 

+ E E + ^?Pil, k;l,k + l-t)+ZEil + ^?P{h k;l + l, k; tjMl) 

k=0 1=0 k=0 1=0 J 
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Then, using the master equations ()A1I - IX6)) . it can be shown that 
d 



- wt) E {Cl, + Cl^,., + Cl,^, + Cl,} - {wD - wt){C% + Clo) 

k=0 

+ 1{u + wt) + liwD - WT){b% + bl,o)] ■ (A42) 



hm ^ < > =d^ 
t^Qo at 



Also, the following equation can be derived using Eq. IA391 

, +00 
^hm - (< l{t) >2) = 2d' A ^ [Cl, + CO+i,, + C,\,+i + Cl,] . (A43) 

k=0 

The formal expression for dispersion is given by 

D = ]-]im^(<l{tf>-<l{t)>'). (A44) 
2 t^oo at 

Then, substituting into this expression Eas. IA421 and IXl3l we obtain 

^ I fc=0 

\{u + wt) + \{wd- WT){bo,o + bi o)} ■ (A45) 

Note, that Tqq = Tqq = Tq and for sum of all Tf^^m we have from Eq. IA38I 

Finally, after some algebraic transformations of Eas. lA37l and lA451 we derive 
the final expression for the growth velocity, V and dispersion, D, which are 
given in Eqs.na'nd[21in Section II. Note that the constant Tq cancels out in 
the final equation. 
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The mean size of ATP-cap can be calculated as 

+00 +00 
< iVcap >= Yl + ^ik) + (2fc + + bl,^,) = kq\l - q) = 

k=0 k=0 

(A47) 

The average relative fluctuation in the size of the ATP-cap, by definition, is 
given 



< Nip >-< N,ar, 



^ cap ^ ^ J V cap ^ 



(A48) 



where 



k=0 

.2 Q„3 



k=0 ^ ^' 

Then from Eg s . I A49l and EZTI we have, 

+ 3(l-g) > 1. (A50) 



a' 1 



< Neap >2 q 
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Figure Legends 

Figure ^ 

Schematic picture of polymer configurations and possible transitions in the 
vectorial model of the single actin filament's growth. The size of the monomer 
subunit is d, while a is a shift between the parallel protofilaments equal to 
a half of the monomer size. Two protofilaments are labeled 1 and 2. The 
transition rates and labels to some of the configurations are explained in the 
text. 

Figure [2J 

Comparison of the growth velocities for the barbed end of the single actin 
filament as a function of free monomeric actin concentration for the random 
and for the vectorial ATP hydrolysis mechanisms. A vertical dashed line 
indicates the transition point c' (in the vectorial hydrolysis) . It separates the 
dynamic regime I (low concentrations) from regime II (high concentrations). 
The kinetic rate constants used for calculations of the velocities are taken 
from the Table 1. 

Figure „ 

Dispersion of the length of the single actin filament as a function of free 
monomer concentration for the barbed end and for the pointed end (the 
vectorial mechanisms). The kinetic rate constants are taken from the Table 
1. Vertical dashed lines indicate the critical concentrations Ccrit for the 
barbed and the pointed ends; thin solid line corresponds to the concentration 
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of the treadmilling, which is also almost the same as the transition point 
for the barbed end. The transition point for the pointed end is at 0.85 /iM. 
Total dispersion is a sum of the independent contributions for each end of 
the filament. 

Figure ^ 

The size of ATP cap as a function of free monomer concentration for the 
barbed end of the single actin filament within the random and the vecto- 
rial mechanisms of ATP hydrolysis. Thick solid lines describe the vectorial 
mechanism, while dotted lines correspond to the random mechanism. The 
kinetic parameters for constructing curves b and d are taken from the Table 
1. For curves a and c the kinetic rate constants are the same with the excep- 
tion of the smaller hydrolysis rate = 0.03 s~^. Vertical dashed line and 
thin solid line indicate the transition point c' and the critical concentration 
Ccrit, respectively, for the curve b. 

Figure 

The fractions of the capped configurations q for the barbed end of the single 
actin filament as a function of free monomer concentration. The results for 
two mechanisms of hydrolysis are presented. Thick solid lines describe the 
vectorial mechanism, while dotted lines correspond to the random mecha- 
nism. The kinetic parameters for constructing curves b and d are taken 
from the Table 1. For curves a and c the kinetic rate constants are the same 
with the exception of the smaller hydrolysis rate = 0.03 s^^. 
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